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Abstract. The separate universe conjecture states that in General Relativity a density per¬ 
turbation behaves locally (i.e. on scales much smaller than the wavelength of the mode) as a 
separate universe with different background density and curvature. We prove this conjecture 
for a spherical compensated tophat density perturbation of arbitrary amplitude and radius in 
ACDM. We then use Conformal Fermi Coordinates to generalize this result to scalar pertur¬ 
bations of arbitrary configuration and scale in a general cosmology with a mixture of fluids, 
but to linear order in perturbations. In this case, the separate universe conjecture holds for 
the isotropic part of the perturbations. The anisotropic part on the other hand is exactly 
captured by a tidal field in the Newtonian form. We show that the separate universe picture 
is restricted to scales larger than the sound horizons of all fluid components. We then derive 
an expression for the locally measured matter bispectrum induced by a long-wavelength mode 
of arbitrary wavelength, a new result which in standard perturbation theory is equivalent to a 
relativistic second-order calculation. We show that nonlinear gravitational dynamics does not 
generate observable contributions that scale like local-type non-Gaussianity f^f, and hence 
does not contribute to a scale-dependent galaxy bias A6 oc k~‘^ on large scales; rather, the 
locally measurable long-short mode coupling assumes a form essentially identical to subhori¬ 
zon perturbation theory results, once the long-mode density perturbation is replaced by the 
synchronous-comoving gauge density perturbation. Apparent contributions arise 

through projection effects on photon propagation, which depend on the specific large-scale 
structure tracer and observable considered, and are in principle distinguishable from the local 
mode coupling induced by gravity. We conclude that any observation of beyond these 
projection effects signals a departure from standard single-clock inflation. 
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1 Introduction and summary of results 

An accurate theoretical description of the nonlinear large-scale structure is one of the ma¬ 
jor goals of theoretical cosmology. Although it is by nature a difficult problem, nonlinear 
structure formation is a key ingredient toward decoding the initial cosmological perturba¬ 
tions in the primordial Universe and testing the large-scale dynamics of gravitation. This is 
especially important given the plethora of current or incoming large-scale structure surveys 
aiming at measurements of unprecedented accuracy. A general relativistic treatment is essen¬ 
tial if perturbations on scales close or comparable to the horizon are being considered. Since 
the initial conditions of structure (e.g., due to inflation) are generally specified when they 
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are outside the horizon, this issue necessarily arises when attempting to connect structure in 
the observable universe to the initial conditions. 

One of the primary reasons of why we would like to connect late-time observables of 
large-scale structure to the initial conditions at nonlinear order is the study of interactions 
during inflation, in particular through primordial non-Gaussianity. Specihcally, the interac¬ 
tion between long-wavelength and short-wavelength perturbations is an interesting diagnostic 
of inflationary physics, as it can distinguish between one or more light degrees of freedom 
during inflation. This can be most directly measured from the squeezed configuration of 
the three-point (bispectrum) or higher-point functions, where “squeezed” means that one 
Fourier mode has a wavelength much longer than all the other modes. For single-clock in¬ 
flation, model-independent squeezed-limit consistency relations are known to hold [1-9], the 
simplest representative of which reads = 5/12(ns — 1), where /^l is the amplitude of 
local-type non-Gaussianity in the primordial (Newtonian gauge) potential. Physically, they 
encode the existence of a single preferred clock. A detection of local-type primordial non- 
Gaussianity that departs from the unique single-field prediction would prove the existence of 
more than one light field during inflation. 

In terms of the matter three-point function, local-type non-Gaussianity schematically 
leads to a squeezed limit of ^ 

(5(ki),5(k2)5(ki))' oc {aH/ksf PgikL) Ps{ks), where k^ = (ki - k 2)/2 , (1.1) 

where corrections are suppressed by {kr/ks)'^ and we have assumed ki < /ceq so that the 
transfer function is unity. While consistency relations are a statement about transformation 
properties under diffeomorphisms, in a local observer’s frame the diffeomorphism freedom 
is completely fixed by the observer’s proper time and distance standards. As shown in [10] 
(see also [11-14]), this implies that squeezed-limit consistency relations drop out of locally 
measurable quantities. Non-single-held models which violate the consistency relation will 
on the other hand contribute to physical observables and can be measured in LSS surveys 
through the bispectrum or scale-dependent bias b{k) oc k~‘^ [15]. 

Thus, an important question to ask is whether nonlinear gravitational evolution of the 
matter density held, neglected in Eq. (1.1), will lead to contributions of lo matter 

statistics. Investigating this requires a relativistic treatment of perturbations at second order 
since both the long- and short-wavelength modes have to be followed through horizon crossing. 
This has been done for example in [16-18]. Unfortunately, standard higher-order perturbative 
approaches for nonlinear mode-coupling are often plagued by unphysical gauge artifacts due 
to the freedom to parameterize perturbations in different coordinate systems (e.g., [19]). 

An alternative approach to this problem is based on the intuition that a local observer 
who has access only to short comoving distances ~ 1/ks, which are much smaller than the 
scale of variation of the long-wavelength modes, would interpret local small-scale physics 
within a Friedmann-Lemaitre-Robertson-Walker (FLRW) background, i.e. a “separate uni¬ 
verse” [20]. To formulate this intuition rigorously, one starts with a coarse-grained universe 
in which small-scale inhomogeneities with ks greater than a fixed comoving scale A are first 
smoothed out. One then studies how small-scale structure with ks > A. evolves in this “local 
background” spacetime, instead of the “global” background universe. To solve small-scale 

prime on the bispectrum indicates that a prefactor of (27r)®(5D(ki -f k 2 -I- ki) imposed by statsitical 
homogeneity has been dropped. We will use this notation in the rest of the paper. For the moment we will 
disregard the issue of which gauge 5(k) is to be evaluated in. 
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clustering in this modified background, one can take a Lagrangian perspective — if non- 
gravitational forces are negligible at large distances, one can follow inertial observers in the 
coarse-grained universe and their vicinity. As initially proposed in [10] and rigorously proven 
in [21], one can locally construct a frame. Conformal Fermi Coordinates (CFC), valid across 
a comoving distance of order at all times, which realizes this physical picture. In CFC, 
the observer interprets the small-scale structure around her as evolving in an FLRW universe 
that is modified due to long-wavelength modes with comoving wavenumber ki < A, with 
corrections quadratic in the spatial distance to the observer 


ds'^ 


o?p{tf) 


—-|-(1-|-Xp/4) -|-0[(xp)^] , 


( 1 . 2 ) 


where ap is the local scale factor and Kp is the local spatial curvature. The fact that the 
metric is valid over a hxed comoving scale is very important. This means that we can apply 
the CFC frame even when the small-scale modes are outside the horizon. This is crucial in 
order to connect to the initial conditions given by inflation. Therefore, within CFC one can 
keep track of any general relativistic effect that could arise as either long- or short-wavelength 
perturbations cross the horizon. 

This “separate universe” approach enjoys at least two advantages. On the one hand, 
gauge freedom is eliminated because the calculation directly yields what the CFC observer 
physically measures. On the other hand, it easily relates to small-scale structure formation 
in homogeneous FLRW backgrounds, which we know very well how to handle. An important 
requirement for this is that ap{Tp) indeed corresponds to a locally observable scale factor. 
This can be achieved easily by using the local velocity divergence as described in [21] and 
Sec. 3 below. 

The separate universe picture strictly holds if the local spacetime is indistinguishable 
from an unperturbed FLRW universe within a small region of space but for an extended (or 
even infinite) duration of time. We find that this is in general only true if the observer and 
all components of matter with equation of state w ^ —1 are comoving. This implies that the 
separate universe picture is restricted to scales larger than the sound horizon of all relevant 
fluids. 

One of the main results of this work is the precise form of the CFC metric Eq. (1.2) and 
the relation of the local scale factor ap and curvature Kp to the long-wavelength perturba¬ 
tions, which we assume to be linear. The latter relation becomes simplest in synchronous- 
comoving (sc) gauge, where we obtain (for a flat global cosmology) the surprisingly concise 
relations 


apiTp) 

a{Tp) 


= I-3A. 


J ^sc(r) 

RTi?) and 


-”o 


(1.3) 


where Age is the long-wavelength density perturbation, a(r) is the scale factor of the back¬ 
ground universe, and D{t) is the linear growth factor normalized to a(T) during matter 
domination. Furthermore, the tidal corrections 0[{xy)‘^] in Eq. (1.2) can be written as 


= aUrF)l - [1 + d4 + [1 - 23-^1 \ (j 

i (l + iFi^xl/d) J 

where ^p, ^p are simply related to the scalar potentials in the conformal Newtonian gauge 
through 

^ (^dkdi<^ - XpXp , (1.5) 
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and an analogous relation holds for 'I'j? in terms of 'h. Note that involve the trace- 

free part of the second spatial derivatives of <1> and i.e. the Newtonian tidal tensor, and 
thns vanish for an isotropic perturbation. 

These results can be summarized in words as: 

On scales larger than the sound horizon of the fluid, the effect of a long-wavelength mode as 
measured locally is completely captured by a modified local scale factor and spatial curvature, 
Eq. (1.3), and a pure tidal field, Eq. (1.5). 

Ref. [10] showed that the consistency relation for single-clock inflation, transformed 
into CFG, states that there is no primordial correlation between long- and short-wavelength 
potential perturbations. This provides the trivial initial conditions for the small-scale fluctu¬ 
ations in CFG. Integrating the growth of small-scale density perturbations 5 in the presence 
of a long-wavelength mode A(ki), we then obtain the leading contribution from nonlinear 
gravitational evolution, valid for arbitrarily small k^, 

5 e = ^Asc + • (1-6) 

Here, <5^^^ is the linear small-scale density field (i.e. in the absence of long-wavelength modes), 
while Asc is the long-wavelength density perturbation (in synchronous gauge), and = 
{didj/d‘^—6ij/3)Asc is proportional to the tidal tensor. Eq. (1.6) shows that the locally observ¬ 
able small-scale density in single-clock inflation has no contributions that scale as local-type 
non-Gaussianity (this holds not just for quadratic f^f but equally for higher order terms 
5 j^£,...). This implies that there is no scale-dependent bias of large-scale strncture tracers 
even when taking into account nonlinear gravitational evolution fully relativistically. We fur¬ 
ther derive the leading squeezed-limit matter bispectrum from Eq. (1.6) [Eq. (5.44)] which 
proves this point. Specifically, this is the bispectrum which would be seen by a “central 
observer” if synchronons observers distributed on the past lightcone of the central observer 
communicated their local densities and power spectra of small-scale fluctuations. This bis¬ 
pectrum is suppressed over local-type non-Gaussianity by a factor of . 

In order to provide predictions for correlation functions measured from Earth, we need to 
inclnde photon propagation (“projection”) effects. These correspond to mapping the locally 
measnrable observables in physical space to the observer’s coordinates of measured photon 
redshift and arrival direction. These effects do contain contributions that correspond to order 
unity f^f. However, we stress that projection effects are purely kinematical, not dynamical, 
and depend on the details of the tracer considered. Thus, they are in principle distinguishable 
from the locally measurable local-type non-Gaussianity, which, as we have shown, is not 
generated by gravitational evolution and must come from the initial conditions. A detailed 
treatment of projection effects is beyond the scope of this paper. Eor recent literature on 
this subject, see [10, 22-24]. 

The remainder of this paper is organized as follows. In Sec. 2, we demonstrate the con¬ 
cept of “separate nniverse” by presenting a fnll general relativistic proof that a compensated 
spherical top-hat region of overdensity embedded in a ELRW universe behaves exactly as a 
separate ELRW universe (with different spatial curvatnre). In Sec. 3 we specialize to the 
GEG for a free-falling observer in the presence of a long-wavelength scalar metric perturba¬ 
tion at linear order. We discuss in Sec. 4 the sufficient and necessary conditions for which 
local ELRW expansion is exact, i.e. for which the spacetime is truly locally indistinguishable 
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from an FLRW universe. The explicit calculation of the locally measurable squeezed mat¬ 
ter bispectrum for an Einstein-de Sitter (EdS) universe {Vtm = 1) with single-clock initial 
conditions is presented in Sec. 5. Concluding remarks are given in Sec. 6. 

Regarding our notation, up to including Sec. 3 we distinguish between the global co¬ 
moving position X* and the CEC comoving position so as to highlight the distinction 
between the global coordinate system and the CEC. After that we will exclusively use spatial 
CEC coordinates and denote them simply by x*. On the other hand, we always explicitly 
distinguish between the CEC times Tp ,tF and the global times r, t. 


2 Proof of the separate universe conjecture for a compensated tophat 


The “separate universe conjecture” states that a spherically symmetric perturbation in an 
ELRW background (taken to be ACDM in this section) behaves like a separate ELRW uni¬ 
verse with different matter density and curvature. This holds up to higher spatial derivative 
corrections, but at all times. We now provide a proof of this statement for a spherical com¬ 
pensated tophat, or just tophat for short (see hgure 1), without making any assumptions 
about the amplitude or wavelength of the perturbation (i.e. we do not assume that it is 
subhorizon). Higher derivative corrections are avoided by assuming a homogeneous density 
perturbation. 

Consider a “background” ELRW universe described in comoving spherical coordinates 

{t,ro,8,(t>) by 

=—dt'^ + a'^{t) [{I — Kr^) + , ( 2 . 1 ) 

where we have allowed for spatial curvature K. The label “o” stands for “outside”, as it will 
become clear shortly. Note that a{t)ro is the area radius, so that the proper surface area of a 
2-sphere of radius Tq is 4:7ra‘^{t)r‘^.‘^ This will become very useful in the following. We assume 
that the universe is filled with dust of uniform density p{t) and a cosmological constant. 
Now consider cutting out a sphere of comoving radius Xo = Ro^c from the ELRW background 
(where “c” stands for “comoving”), and collapsing the matter in the sphere into a point mass 
M. If we choose Ro,c to be constant in time, then M is conserved.^ The authors of [25] 
first showed for a matter-dominated universe that there is a unique solution consisting of a 
Schwarzschild metric interior to Ro,c smoothly matched to the ELRW background Eq. (2.1) 
at Xo = Ro,c for all times t (more precisely, the metric is continuously differentiable at the 
boundary). This result has since been generalized to include a cosmological constant [26, 27], 
in which case the point mass metric becomes the Schwarzschild-de Sitter solution: 


dssds — “ 
V{rs) = 1 


V{rs)dtl -t- V ^(rs)dr^ + r^dfi^ 

2GM 1 . 9 
-xAr: , 


( 2 . 2 ) 


where the label “s” stands for “Schwarzschild”. Hence, the matter surrounding the shell 
is oblivious to the fact that the interior has been collapsed to a black hole, while a test 
particle anywhere inside Ro,c is oblivious to the surrounding homogeneous matter. Note 

^The geodesic radius on the other hand, i.e. the proper radius along a radial geodesic, is given by y, where 
X is defined through Xq = s\nK{^\K\x)/\J\K\ and sin^ x = sini, x, sinhx for If > 0, = 0, < 0, respectively. 

^In a fiat background, this mass is simply (A-k/Z) while for K ^ Q, the mass-radius relation 

is more complicated. Independently of curvature, the interior mass scales as a^{t) for constant i?o,c and is 
conserved. 
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Figure 1: Illustration of the setup used to prove the separate universe approach. 

that by matching the angular part of the metric, we can immediately identify = a{t)ro. 
Moreover, is the area radius for ^ 53 ^ 3 . In the interior coordinates, the boundary is thus 
set by Ts^bound = Ro{t)- Using this relation, it is then easy to show that at the boundary 
the geodesic equation derived from the metric Eq. (2.2) exactly matches that given by the 
FLRW background Eq. (2.1), which is simply Ro{t)/a{t) = Rq^c = const, once the Friedmann 
equation for a{t) is inserted. This is of course a consequence of the metric being continuously 
differentiable at the boundary. 

Consider now the case where we do not collapse the matter inside Rq into a black hole, 
but compress it to a finite radius = Ri{t) (where “z” stands for “inside”), maintaining a 
homogeneous density p (see Fig. 1). Birkhoff’s theorem [28], generalized to include a cosmo¬ 
logical constant, states that the unique spherically symmetric vacuum solution to Einstein’s 
equation is Schwarzschild-de Sitter. Thus, Eq. (2.2) still describes the spacetime outside of 
the mass, Vg > Ri- We now perform the exact same matching as derived by [25-27], but 
inverting outside and inside. In fact, nothing in the matching is particular to the case of vac¬ 
uum inside a homogeneous matter distribution; it equally applies to a homogeneous density 
distribution inside vacuum. Thus, there exists a unique FLRW solution of the form Eq. (2.1) 
with density p that smoothly matches to the Schwarzschild-de Sitter metric at Ri, where p 
is determined by mass conservation. 

Moreover, we can use the geodesic motion that follows from the metric Eq. (2.2) to derive 
the scale factor a{t) in the interior FLRW solution. The geodesic equation for Eq. (2.2) and 
purely radial motion, = dxg/dt = {vP,u'') where t is the proper time, becomes 


dvP ldV{rs{t)) 
dt 2 dt 


(2.3) 


which can be integrated to give 


Ri{t) -|- V{Ri{t)) = const, 

where we used Ri = dRi/dt = vP. Inserting Eq. (2.2), we obtain 


(2.4) 



2GM 1 


-I— KR^ + const. 
3 


(2.5) 


R^ 












This can be turned into a Friedmann equation by introducing a scale factor a{t) oc Ri{t) 
different from the “background” a{t). The normalization of the scale factor is arbitrary, and 
we can choose to normalize to M = {4:7r/3)p{t)Rf{tQ)a^{t) where to is some reference time 
and p{t) oc Dividing Eq. (2.5) by i??, we obtain 


= 



8ttG _, , 1 , 

-I^pM + 3A 



(2.6) 


We have relabeled the integration constant suggestively as —Ki. This is indeed the Fried¬ 
mann equation for a ACDM universe with background matter density p and curvature K^. 
Assuming that the perturbation was initialized with a very small amplitude at an early time 
tj when the effect of A is negligible, one can then subtract the background FLRW equation 
following from Eq. (2.1), and linearize in {Ri — Ro)/Ro at U. As shown in Sec. 4 below (see 
also [29]), this yields a relation between the initial (linear) density perturbation and curvature 

K, = , (2.7) 

where for simplicity we have assumed that the background ELRW universe is flat (although 
this is not a necessary assumption). 

Thus, we have proved that, at least for the case of a compensated tophat density profile, 
a spherically symmetric perturbation in a ACDM background evolves exactly as a separate 
curved ACDM universe. We have not assumed that the perturbation is small or that the 
scales Ro, Ri are subhorizon. Indeed, neglecting A for simplicity, the Schwarzschild-de Sitter 
metric in the vacuum surrounding the perturbation cannot be perturbatively approximated 
as Minkowski if 

~ 1, hence {HiRif ~ 1, (2.8) 

rii 

since Gp ~ Hf. Thus, for a horizon-scale perturbation, the vacuum exterior is far from 
Minkowski. Nevertheless, the separate universe description is exact. In fact, an observer in 
the vacuum region surrounding the overdensity would see this separate universe as a black 
hole. 

One might wonder how this can be applied to an underdensity (void) rather than an 
overdensity. In this case, the inner FLRW solution which we have called p now becomes 
the background p, while the outer, less dense FLRW solution becomes the “perturbation” p. 
From the perspective of an observer inside either FLRW solution, the situation is completely 
symmetric, as they cannot tell by any local measurement whether they are embedded in a 
larger “background” universe. 

In Sec. 4, we will return to the separate universe picture in the context of general scalar 
perturbations in a cosmology with multiple fluids. 


3 Scalar perturbations in the CFC frame 

In this section, we describe how long-wavelength adiabatic scalar perturbations are treated 
in the CFC frame. This provides the basis for reconsidering the separate universe ansatz in 
the following sections. 


- 7 - 





3.1 Review of Conformal Fermi Coordinates (CFC) 

We first briefly review the basic concept of Conformal Fermi Coordinates (CFC), which was 
first introduced in [10] and defined rigorously in [21]. Consider an observer free-falling in some 
spacetime. We are mostly considering applications in the cosmological context, in which this 
spacetime is approximately FLRW, but this does not have to be the case. Her trajectory is 
a timelike geodesic G, whose tangent vector we will denote as eo- One can construct a frame 
where the spatial origin is always located on this geodesic, and in which the metric takes the 
form^ 

x^p)] , = 0[{x^pf] . (3.1) 

Thus, in these coordinates the metric looks like an unperturbed FLRW metric in the vicinity 
of the observer’s trajectory at all times, up to corrections that go as the spatial distance from 
the origin squared. Note that a metric of the form Eq. (3.1) means that an observer at the 
origin x'^p = 0 is free-falling, with proper time tp given by 

tp{Tp) = [ ap{T')dT'. (3.2) 

Jo 

Briefly, this frame can be constructed as follows. For now, let ap{x) be a positive scalar field 
in the neighborhood of G. We introduce the conformal metric 

giiAx) = aF^{x)gf,u{x ). (3.3) 

We then define a “conformal proper time” Tp along G through 

dtp = ap{P)dTp , (3-4) 


where P denotes the point on the central geodesic which has proper time tp. This defines 
our time coordinate = Tp. At some point P along G we can construct an orthonormal 
tetrad (ea), where cq = d/dtp is again the tangent vector to G, with {ea)^{ebYg^y = rjab. 
The tetrad is defined at all other points on G by parallel transport, so that this condition 
is preserved. A point Q corresponding to CFC coordinates Xp is then located as follows. 
First, we move to the point P on G specified by x^. We then construct a spatial geodesic 
h{xp-, x^p; A) of the conformal metric which satisfies h{x^‘, Xp] 0) = P, and whose tangent 
vector at P is given by 


^h^{xY,xYA) 


A=o 



s = {SijX^pX^pY/'^ . 


(3.5) 


The point Q is then found by following this geodesic for a proper distance A = ap{P)s. As 
we show in [21], this yields a metric in the form Eq. (3.1). Moreover, the leading corrections 
are given by 


^00 — ~ 

- -- 

hfj= - ^RujmXpXf , (3.6) 

^Here we have let qf be a function of the CFC time coordinate tf only. While this is not strictly necessary 
(see [21]), this is the most natural choice as constant-scale-factor surfaces then coincide with constant-proper- 
time surfaces. 


Roiim^F^^ 
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where is the Riemann tensor of the conformal metric, evaluated in the CFG frame at 

x^p on the central geodesic. 

Note that if we now were to choose ap = 1, the CFG would reduce exactly to the 
ordinary Fermi Normal Goordinates [30]. Then, the metric perturbations would contain 
the leading corrections due to Hubble flow oc restricting the validity of the coordinates 

to subhorizon scales. Instead, we will choose ap so that it captures the locally observable 
expansion of the Universe, which automatically extends the validity of the coordinates to the 
spatial scale of perturbations which can be superhorizon. The well-defined prescription that 
fixes ap up to a multiplicative constant will be given in Sec. 3.2.1. 

3.2 Scalar perturbations around an FLRW universe 

We now turn to an FLRW universe with scalar perturbations in the conformal-Newtonian 
(cN) gauge. 


ds^ = o?{t) [—(1 -I- 2<l>)dT^ + (1 “ 2^)5ijdx^dx^]^ . (3.7) 

In the absence of anisotropic stress we have T = <I>; however, we will not make that assump¬ 
tion in this section. Before we apply the GFG construction, we note that we are implicitly 
performing a coarse-graining of the metric Eq. (3.7) on some scale A. This is because in the 
actual universe, scalar perturbations exist on all scales, so that without any coarse-graining 
the metric Eq. (3.1), which assumes that the corrections are small, is only valid on in¬ 
finitesimally small scales. In the following, <I> and 'L will thus denote coarse-grained metric 
perturbations. 

Further, we will work to linear order in <I> and 'k. The perturbative expansion in <k, 'k 
should not be confused with the power expansion in The former expansion is valid as 
long as |<k, T| <C 1, and is chosen here for simplicity; the GFG construction also works for 
spacetimes that differ strongly from FLRW. The latter expansion on the other hand is good if 
is smaller than the typical variation scale for 'k, which is the fundamental expansion 
parameter in GFG. Following standard convention, spacetime indices are assumed to 

be raised and lowered with the metric while latin indices i,j, k ,... are raised and lowered 
with Sij. Gorrespondingly, denotes the flat-space Laplacian 

92 = d^^didj . (3.8) 


3.2.1 CFG construction 

We begin with deriving the tetrad. Gonsider a free-falling observer traveling along the time¬ 
like central geodesic G. His 4-velocity is given by 

= (eo)^ = (1 - r) , (3.9) 

where the 3-velocity U* is considered as first-order perturbation.^ U* obeys the geodesic 
equation, given by 


w' + nv^ = , 


(3.10) 


^Throughout, L* will denote the large-scale peculiar velocity in global coordinates, while we will later 
encounter the small-scale peculiar velocity u*. 
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where a prime denotes derivative with respect to conformal time r. Our choice of the spatial 
components of the orthonormal tetrad is 

(ejf = + , (3.11) 

which we have aligned with the global coordinate axes for simplicity and without loss of 
generality. 

We now want to determine ap hy defining the locally observable expansion rate of 
the Universe to match the CFG Hubble rate Tip ^ ap^dap/drp. Following [21], consider 
the convergence of the geodesic congruence along the central geodesic G. This is a 

locally measurable quantity, which intuitively corresponds to the change in time of the proper 
volume of a bundle of geodesic trajectories in the neighborhood of the central geodesic G. In 
an unperturbed FLRW universe, we have 


= 3H = 371/0,. 


(3.12) 


Thus, we define at all points along the geodesic G 




G 


(3.13) 


In terms of cN gauge quantities the convergence of a geodesic congruence is given by 

= 3H[l-^- + {37i)-^diV^] . (3.14) 


The local Hubble parameter is then given by 


Hp = H 






(3.15) 


Note that Ref. [35] in their construction using the standard Fermi Normal Coordinates (FNC) 
only include the last, velocity-divergence term, which dominates when the long-wavelength 
mode is deep inside the horizon. However, the second term here is a necessary correction 
on (super-)horizon scales, since on those scales W does not correspond to physical peculiar 
velocities. 

Consider, for example, the motion of two nearby fluid elements in the presence of a 
superhorizon long-wavelength mode. For adiabatic initial conditions, the two fluid elements 
will then have vanishing relative peculiar velocity, which scales as (k/T-l)'^^, i.e. the fluid 
elements have a constant, infinitesimal coordinate separation Ax. The physical separation is 
dphys = a(l — T) Ax. The proper time interval is a(l -|- <h)(ir. The physical relative velocity 
Uphys, i-e. the rate of change of the proper distance between them, then reads 


'l^phys 


I 


d 


a(l -|- 4') dr 


^phys — 


1 


d 


a(I -|- 4)) dr 


[a(l - T) Ax] = ?^ (1 - $ - T - —T' 


7i 


= H 




o(l - T) Ax = Hpdphys 


Ax 

(3.16) 


where the last equality holds since we neglect terms of order (fe/'R)^<h. Therefore, the relative 
velocity Uphys and Hp satisfy exactly the local version of Hubble’s law, showing that Eq. (3.15) 
is indeed the proper expression for the local Hubble rate. 
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Since so far we have only fixed Hp = dlnap/dtp, the local scale factor ap is defined 
up to a multiplicative constant. In other words, we have a residual freedom to rescale 
{ap, Tp} —)■ {oirA, Tp/X}, which keeps the proper time tp unchanged. Of course, this constant 
is arbitrary and should cancel out of any proper observable. We will here fix the constant 
in order to make our results of the next sections more transparent. Specifically, we demand 
that 

... 1 as Ti;’ —)■ 0, (3. 

a[Tp) 

which means that at early times, the local scale factor-proper time relation is the same as 
that in the unperturbed background cosmology. The ratio of scale factors at a fixed spacetime 
point is then at early times given by 



apijp) 
Inn —-—^ 

tf-S-O a[T[Tp)) 


= 1 + 


(3.18) 


where <hini is the asymptote as rj? —)• 0 of the potential along the geodesic. A direct integration 
gives 


^ = 1 + + r dr' (-<!>' + v) , (3.19) 

a(r) 3 Jo V 3 y 

In order to obtain the leading corrections to the CFC metric (either via Eq. (3.6) or an 
explicit coordinate transformation of the metric), we also need the second derivative a'p. 
This can be obtained easily from Eq. (3.13), 


1 dTipi^Tp) 1'Hp{tp) 

“I’('^f) drp 3 ap{Tp) 

where again all quantities are evaluated on G. 


V;.t/^ + ^(eo)"VFV^[/^ 


(3.20) 


3.2.2 CFC metric 

We can now derive the remaining corrections to the CFC metric, making use of the 
geodesic equation for F*. This yields in terms of the conformal Newtonian gauge perturba¬ 
tions 

/loo = - (dkdi^ - XpXp (3.21) 

doi — 2 {Td ~ {dkiyi ~ diiVk) XpX^p - {dkidi — dkidi) (T -|- XpX^p (3.22) 
hfj = ^ {6ijdidk'i> + diidjdk'it - dijdkdi'it - dkAdj'it) XpXp (3.23) 

{6ij6ki — dikSji) XpX^p . 

y 

Let us consider a spherically symmetric configuration about the CFC observer. It amounts 
to setting W, di^ and to zero, and replacing didj^ with {l/3)dijd‘^^ and the same for 
T. We then find that /i^ and /i^ vanish, while hfj reduces to 

hfj = ^Kf {xp,ixpj - 6ij rp) , (3.24) 
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where we have suggestively introduced 


Kf=‘^ - ndiV^) . (3.25) 

O 

The tensorial structure oc (^XpX^p — distorts the proper distance only along directions 

perpendicular to the radial direction, as guaranteed by our construction of CFC (recall that 
rp is defined as the geodesic distance in the conformal metric). We are free to choose a 
different radial coordinate (see App. A for a discussion of the residual gauge freedom in the 
spatial component of the CFC metric) 

rp —^1 - , drp —^ “ \^F drp . (3.26) 


In terms of this new radial coordinate, the area radius of the spacetime (see Sec. 2) is given by 
{1 + Kp rp/4:)~^rp. Under this rescaling, the CFC metric simply becomes (valid to 0[(x);.)^]) 


ds'^ 


ap{Tp) 



dijdxpdx^p 

(l+Kpr^p/A)^ 


(3.27) 


The spatial part of the metric is the familiar stereographic parameterization of a curved, 
homogeneous space. 

We now go back to general anisotropic configurations of T. In the next section, 
we will see that, in the context of scalar perturbations considered here, it is necessary that 
all cosmic fluids are comoving with velocity V) for the separate universe picture to hold. 
The comoving condition implies that in fact vanishes, by way of the time-space Einstein 
equation. Therefore, we will set to zero in the remainder of the paper. We can again 
use the residual gauge freedom described in App. A.2 to bring hC into a more familiar form. 
After some algebra, this yields 





xWf 

XpXp 


drp 

5ijdx''pdx^p I 

{1 + Kpr‘^/Af> 


(3.28) 


We have thus put the CFC metric into the conformal Newtonian form. The metric per¬ 
turbations <I>, T enter in two distinct ways. First, the isotropic part of the perturbation 
(proportional to and <I>, T themselves as well as their time derivatives) modifies the 
background apijp) and leads to spatial curvature. The anisotropic part, which is completely 
determined by 

(3.29) 

and the analogous Kf^, enters as a tidal field. Note that since the CFC frame was constructed 
based on local observables, no gauge modes yield any observable imprint in the CFC metric. 
Further, it is important to emphasize that no subhorizon assumption has been made about 
either <I>, T or the CFC metric itself; Eq. (3.28) is valid on arbitrarily large scales as long as 
the corrections are small. We will turn to the interpretation of the CEC metric in the next 
section. 
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4 Separate universe revisited 


In the previous section we derived the CFG frame metric for scalar perturbations, which 
appears to be of the FLRW form with tidal corrections. However, in order to test whether an 
observer along the central geodesic would indeed interpret local measurements as an FLRW 
universe in the absence of tidal corrections, we need to derive the equation for ap and compare 
it to the Friedmann equations. Further, we need to verify that the “curvature” Kp is indeed 
constant. Throughout we will set the tidal corrections in Eq. (3.28) to zero in keeping with 
our linear order treatment. 

Let us consider an observer traveling along the central geodesic who performs local 
measurements. Suppose further this observer knows about General Relativity. She will 
compare the local Hubble rate Hp with the Friedmann equations, the first of which reads 



SttG 



(4.1) 


where pp \s the local rest-frame matter density, coarse-grained over the region considered as 
GFG patch. Using Eq. (3.15), Kf is given in terms of global quantities by 

Kf = + 2% (T' + n^) - hi diV ^, (4.2) 

O 

where Acn is the density perturbation in cN gauge. On the other hand, the spatial curvature 
Kp which can be inferred from, e.g. shining light rays and measuring angles between them, 
is given in terms of global quantities by Eq. (3.25). We thus need to verify whether = 
Kp = const, and under which conditions. Only in that case would the observer actually 
interpret the results of local measurements as an ELRW universe. 

Let us assume that the Universe is filled with multiple uncoupled fluids labeled by 
I = 1,2, •••. Each fluid has a homogeneous equation of state Pi = wjpi, with possibly 
time-dependent wp The homogenous part of each fluid evolves as p'j = —34^(1 -|- wppi, 
which implies 

{n^niy =-{i + 3 wi)n^ni, ( 4 . 3 ) 

where 0/ = pi/p and p is the total (homogeneous) energy density. The acceleration is given 
by 

n'= -\n^{^{i + 2,wi)n^. (4.4) 


In the global cN-gauge coordinates, we let each component have a fractional density pertur¬ 
bation A/ and peculiar velocity V} = d^Vp where Vi is the velocity potential. Note that U* 
denotes the velocity of the central geodesic throughout. We now write down the time-time 
and time-space Einstein equations, which have the sum of all fluids as sources. 




+ = AuhY,^,{i + wi)Vi\. 


(4.5) 

(4.6) 
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Using that 


(4.7) 


Acn = ^ ri/A/, 

I 

Eqs. (4.5)-(4.6) can be inserted into the first two terms of Eq. (4.2) to yield 

Kf = I ( 52 ^ - nd^V^) = Kf . (4.8) 

o 

This is Eq. (3.25), showing that the quantity appearing in the local Friedmann equation is 
indeed the spatial curvature. Thus, Eq. (4.2) provides an equivalent relation between the long- 
wavelength adiabatic perturbation in conformal-Newtonian gauge and the local curvature. 


4.1 Conservation of curvature 

In the FLRW solution, the spatial curvature Kp is constant. This is thus a necessary con¬ 
dition that needs to be satisfied for the separate universe picture to hold. To study this, we 
allow each fluid to have a pressure perturbation 5Vi in addition to the density perturbation 
5pi = piAj. In the rest frame of the fluid, we have 


= cl,I 


(4.9) 


where the rest-frame speed of sound j 


does not necessarily equal the adiabatic sound speed 


eld,I 


P'l 

— =wi - 
Pi 


w, 


3n{i + wi)' 


(4.10) 


because the pressure perturbation can be non-adiabatic. Each fluid then evolves according 
to [31] 


A'j + (1 + wi) d'^Vi - 3 (1 + wi) T' = 3n {wi - cI j) [Aj - 3?7(1 + wi)Vi] 

+3nw'jVi 

c 2 

V/ + (1 - Sell) nvi + —+ $ = 0. 

1 + Wj 


(4.11) 

(4.12) 


Using Eqs. (4.3)-(4.4), Eqs. (4.5)-(4.6), and Eqs. (4.11)-(4.12) (keeping j general), we 
can directly compute the time derivative of Eq. (4.2).® We also assume that U* follows the 
peculiar motion of the non-relativistic matter and therefore satisfies the geodesic equation 
(the cl j = 0 version of Eq. (4.12)). We obtain 


dKp 

dr 


TP f ^iAA - 3%^ f 0/(1 + wi)Vi\ 






^ (1 + wi) nidi [Vi - w) 


L I 


+ (4> - T). 

O 


(4.13) 


The final result is obtained after lengthy but straightforward algebra. The second term 
vanishes if long-wavelength perturbations do not source anisotropic stress <1> = T. 

®Note that d/dr = dr + V'di is in principle the total derivative along the central geodesic; however, the 
spatial derivative term has been neglected as it is higher order in pertnrbations. 
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If the Universe is only composed of non-relativistic matter {w = = 0) plus a possi¬ 

ble cosmological constant {w = —1), and the CFG observer co-moves with the matter fluid 
y = along the central geodesic, Kp is conserved and the observer cannot distinguish 
the spacetime from a curved ACDM universe using any local measurements (that is, mea¬ 
surements over a scale much smaller than the wavelength of the perturbations <h, This 
applies to any isotropic configuration around the geodesic, and is not restricted to the tophat 
considered in Sec. 2. 

If there are extra fluids with other values of wj in the Universe, as postulated in many 
of the alternative cosmologies, (IKp/dr in general does not evaluate to zero. If all fluids 
CO- move with the CFG observer Vj = U*, however, Kp does remain a constant. Since this 
situation requires all fluids to follow geodesic motion, non-gravitational forces have to be 
negligible. This will be true when the sound horizon is much smaller than the wavelength of 
interest, i.e. f a(T)cs(T)dT a/kp [32]. This further generalizes our considerations of Sec. 2. 

Assuming that the GFG observer co-moves with the average cosmic fluid such that 
Vi = 5T^i/{p + V), we can re-write Eq. (4.2) using Eqs. (4.5)-(4.6) to hnd 


Kp 



nd^5T^i\ 
P + V ) 



(4.14) 


where TZ is the gauge-invariant curvature perturbation on comoving slices. Notice that we 
have made no assumption about the decaying adiabatic mode (the second adiabatic solution, 
always present in the small k limit [33]). Its contribution cancels out of this equation. The 
conservation ol Kp is hence directly related to the conservation of TZ and it is valid also 
in the presence of a decaying adiabatic mode. Note that the conditions for TZ = const 
are different than for C, = const, where ^ is the curvature perturbation on uniform density 
slices. The latter is conserved if pressure is only a function of energy density [34], while 
the conservation oi Kp does not exclude non-adiabatic pressure. Thus, TZ, not C; is the 
locally observable curvature, and the separate universe picture relies on conservation of the 
former, not the latter. Note that in [29, 35], the comoving curvature perturbation was 
denoted with (, specifically T^here = — Cthere- Thus, our result Eq. (4.14) agrees with those 
references. Eq. (4.14), together with the discussion after Eq. (3.15) above, thus confirms 
that superhorizon perturbations cannot affect the cosmology as observationally determined 
by measurements within the horizon [36]. 

Finally, one can also verify that the second Friedmann equation is satisfied hy Hp, 
which is shown in App. G. Unlike the first Friedmann equation, this does not impose any 
conditions on the fluid components. Indeed, it is a simple consequence of the Raychaudhuri 
equation applied to the GFG metric Eq. (3.28). 


4.2 Relation to synchronous-comoving gauge 

It is illuminating to connect the curvature Kp, Eq. (4.14) to the synchronous-comoving (sc) 
gauge. Following [38], we write the metric in global coordinates as 


ds'^ = a^(r) 


1 , 


-dr + { -h6ij VVijid 0 


dx^dx^ 


(4.15) 


where TTij = didj — 6ijd'^/3. Here and in the following, we drop the subscript F on the 
spatial GFG coordinates, since the global coordinates will not appear explicitly anymore and 
so there is no ambiguity. The time-time component of the Einstein equation is 


- = dvrGpa^As 

6 2 


(4.16) 
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where Age is the density perturbation in synchronous-comoving gauge. Further, the comoving 
condition y* = 0 implies that = h'. Note that the quantity on the l.h.s. of Eq. (4.16) is 
524 -, i.e. we have 

5^4' = dvrGpa^Asc . (4.17) 


That is, the spatial metric potential in cN gauge is related to the density perturbation in sc 
gauge by the Newtonian Possion equation. This relation holds on all scales, a fact which is 
being used in N-body simulations [37]. 

One can easily verify that the local Hubble rate Hp for the metric Eq. (4.15) is given 
by 

= H+ ]^a-^h'. (4.18) 

Let us now consider a spherically symmetric perturbation in sc gauge. Straightforward eval¬ 
uation of the spatial part of the conformal Riemann tensor, together with didjh = 5ijd'^h(3 
then yields 

h(j = - 0 {xiXj - Sijr^) , (4.19) 

and thus, matching to the spatial curvature (Eq. (3.24)), 


Kp = -^d^{h-0- 

Let us verify again the local Friedmann equation. 


Hl = 


8ttG 


-PF 


Kf 


(4.20) 


(4.21) 


Multiplying the Friedmann equation by a^, using Eq. (4.18), and subtracting the background 
contribution, we obtain 

^Uh' = + ^d^(h - 0 , (4.22) 

which is exactly the time-time component of the Einstein equation Eq. (4.16) (multiplied by 
2/3). 

Let us now restrict to an EdS universe. The continuity equation for pressureless matter 
in sc gauge reads [38] 

A', = -h'. (4.23) 

We thus have 

A'f = —j^a^pAsc — + -^T-iA'^c 

= ^-Hla-^A,, , (4.24) 

where we have used that Age oc a(r) while 'h?' = HQa~^ for EdS universe. Note that for 
any (spatially flat) cosmology where Kp is a constant and which has an epoch of matter 
domination, during which = {Qrno/^m)HQa~^ ^ Q,rnoHQa~^, this expression can be 

simply evaluated during that epoch, so that for any later expansion history (ACDM or other) 


A'f _ 5^ Agc(r) 

H2 - 


(4.25) 
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where D{t) is the linear growth factor normalized to a(T) during matter domination. Thus, 
the local curvature is directly related (during matter domination) to the matter density 
perturbation in synchronous-comoving gauge. Moreover, we find the same relation as in the 
nonlinear compensated tophat of Sec. 2, Eq. (2.7). 


5 Application: the squeezed-limit matter bispectrum on large scales 


Let us now turn to an application of the CFG frame for long-wavelength scalar perturbations, 
namely calculating the physical influence of these perturbations on the growth of small- 
scale structure. We will restrict to non-relativistic pressure-free matter on short scales. As 
described in detail in [21], we coarse-grain the metric on some spatial scale L = A“^, and 
perform the CFG construction with respect to this coarse-grained metric. In [21], we used 
the notation for the long-wavelength part of the metric perturbation, which contains only 
Fourier modes with < A, while the remaining short-wavelength metric was denoted h^, 
with support for > A in Fourier space. Here, for clarity of notation and since we have 
adopted the conformal-Newtonian gauge, we will let T stand for the long-wavelength 
perturbation, while 0, ijj denote the small-scale metric potentials. In the GFG frame, the 
metric is then a simple extension of Eq. (3.28), 


ds^ = a\{TF)l - 


l + 2(j)+ (^kdi^ - x^x'‘ 


dTp 


+ 


1-211)- [ dkdi^ - ) x^J- 


dijdx^dx^ 

{l + Kpr'^/Ay 


(5.1) 


Note that 4>, if) are defined with respect to ap and are thus distinguished from the usual 
parameterisation in global coordinates. The crucial point is however that the GFG frame is 
directly related to observations (see Sec. 6 of [21]), so that this form of the metric is actually 
a desirable feature. For the time being, (p and ip are allowed to be fully nonlinear. 

We now want to solve for the evolution of p, ip in the presence of the long-wavelength 
mode T. At first, one would expect to have to solve the nonlinear Einstein eguations, that 
is, at second and higher order in cosmological perturbation theory. However, as described 
more generally in Sec. 5 of [21], one of the key virtues of the GFG framework is that this is 
actually not necessary, if one restricts to the regime in which small-scale perturbations have 
long entered the horizon (while the wavelength of the long mode is not restricted). Note 
that for adiabatic perturbations, there are no interesting dynamics at all until the short- 
wavelength modes enter the horizon, and the long-wavelength mode, being far outside the 
horizon at this point, cannot have any dynamical impact. 

The reasoning proceeds as follows. First of all, ap and Kp can simply be incorporated 
into the background part of the Einstein tensor. Further, for ks Td we can identify ip = p, 
and only use the scalar Einstein equations, i.e. the 00-component as well as the divergence 
of the Of-component. Finally, we only need to solve the Einstein equations on the central 
geodesic G, i.e. at x* = 0, since we are free to choose the desired fluid trajectory of interest as 
central geodesic, and since we do not need to take additional spatial derivatives of the Einstein 
equations. Glearly, only terms in the nonlinear Einstein tensor with two spatial derivatives 
acting on the or T pieces can contribute to the long-short coupling contribution to the 
Einstein equations on G. Since the Einstein tensor is linear in second derivatives, the short- 
wavelength perturbation in these terms can only contain p, ip without any derivatives (these 
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terms come from contracting indices with the perturbed metric). In other words, the only 
second-order terms in the Einstein equations would have to be of the form 

F{4>) , (5.2) 

and correspondingly with 'h —>■ <h, '0. However, in order for these terms to contribute 

to scalar Einstein equations, the only tensor available to contract with is 6ij. Hence we see 
that these terms have to vanish. 

We emphasize that even the case ks < F. is computable in CFC (as demonstrated in 
Sec. 7 of [21] for the case of long-wavelength tensor perturbations), although the calculation 
becomes technically involved because the trace-free part of the ij-Einstein equation must be 
used and terms of second order in perturbations start to appear in the Einstein tensor. Given 
the limited phenomenological significance of the regime ks <71, we will focus on ks Ti in 
the remainder of this section. This implies that we can set cf) = il) even at second order. 

To summarize: the isotropic part of long-wavelength metric perturbations is entirely 
absorbed into the background apiTp) and Kp, and thus does not need to be dealt with 
explicitly; while the anisotropic part does not enter the Einstein equations for symmetry 
reasons. Thus, in order to solve for the long-short mode coupling in the CFC frame, it is 
sufficient to solve the purely small-scale Einstein equations with a modified scale factor. This 
holds as long as one can restrict to subhorizon short modes which are being influenced by 
long modes of arbitrary scale. The same was found to be true in [21] for long-wavelength 
tensor instead of scalar modes. 

In order to close the system, we need to also consider the fluid equations 

= 0 . (5.3) 

This is again very simple, since we merely have to include a modified background and a tidal 
force which preserves its “Newtonian” form on all scales [Eq. (5.1)]. Again, this is closely 
related to the analogous problem of long-wavelength tensor modes studied in [21]. 

5.1 Outlirre of the calculation 

Let us outline the calculation that will take up the remainder of Sec. 5. We will consider 
the leading coupling of long- and short-wavelength modes at second order in perturbation 
theory, so that we arrive at a second order small-scale solution 

(kp + %) ~ (%), (5.4) 

where indicates linear order solutions. We emphasize that this is merely for simplicity. The 
CFC construction works for both nonlinear large-scale and small-scale metric perturbations, 
and the simplification described above is not restricted to the leading order only. This is one 
of the key advantages of this approach. A higher order calculation in terms of the small-scale 
fluctuations has recently been performed using N-body (“separate universe”) simulations 
[29, 39, 44]. Going to higher order in the long-wavelength fluctuations is beyond the scope 
of this paper and will be pursued elsewhere. As discussed in detail in Sec. 5 of [21], it is then 
sufficient to solve the linear Einstein equations for the long mode, and insert the solution for 
the long-mode metric perturbation into the small-scale fluid equations. 

We will begin by calculating the evolution of small-scale modes in the presence of long- 
wavelength perturbations, i.e. the long-short mode coupling measurable by a CFC observer, 
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separately for the cases of isotropic and anisotropic long-wavelength perturbations (Sec. 5.3- 
Sec. 5.5). We then consider the mapping of CFC coordinates between neighboring observers 
(Sec. 5.6). This leads to a final expression for the locally measurable squeezed matter bis¬ 
pectrum in Sec. 5.7, to be understood in a sense described there. 

The results of Sec. 5.3 to Sec. 5.7 will be derived for a flat Einstein-de Sitter (EdS) 
background. This is merely for calculational simplicity. By replacing the EdS scale factor 
a{t) with the linear growth factor D{t), these results are numerically highly accurate even 
for other cosmologies, e.g. ACDM. 

We conclude this outline by presenting the relevant equations. Since the long-wavelength 
density perturbation is entirely absorbed into the local background apiTp), the only long¬ 
mode contribution to be included in the fluid equations is the local velocity field induced 
by the long-wavelength mode, which we denote as v\^. In CEC [Eq. (5.1)], the equation of 
motion for the latter becomes 

+ '^F + KfjX^ = 0. (5.5) 


This can be integrated to give 




aFiT')Kfj(T')dT' = -- 


apir) Jo 

where we have defined the long-wavelength tidal tensor along G as 


(5.6) 


= 


didj - 




(5.7) 


and the second equality is specific to an EdS universe. 

The 00-Einstein equation, continuity equation, and Euler equation for the small-scale 
modes then become 


5' + di[v^ + {I + 5)vi] = 0 
v[ +%FVi + {VL ■ d)vi + {v ■ d)vL,i = -di(j) , 


(5.8) 


where all spatial derivatives are with respect to x* = Xp. Here and throughout, u* denotes 
the small-scale velocity field, i.e. the total velocity in CEC without the tidal contribution 
Vp due to the long-wavelength mode. Note that we have dropped terms purely quadratic 
in small-scale perturbations in Eq. (5.8) in keeping with our second-order treatment in this 
paper. 

Eq. (5.8) are almost identical to the well known system of equations used to solve 
the evolution of large-scale structure on subhorizon scales. The only impact of the long- 
wavelength mode, apart from the modified cosmology entering via Hf, is by adding the 
tidal-field induced velocity Vp. Note that we have to keep this velocity to order x* in Eq. (5.8) 
since the evolution of the density field on G depends on derivatives of Vp on G. Eor an EdS 
background, the local matter density parameter in CEC is given by 


- 1 = 

™ r,2 u2 

OipFi F 


Kf Kf 


K 


(5.9) 
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where Kp is the local curvature parameter as discussed in Sec. 4. In Eq. (5.8) we have 
neglected terms of the form Kpcf) and replaced the curved-space Laplacian j with the 

flat-space Laplacian. Kp was of order T-lp, these terms would be suppressed on subhorizon 
scales with respect to the remaining terms by the same order of magnitude as terms such as 
'Hpfj), T-Lpct)' which we have neglected. In fact, Kp/'Hjp ~ Age 'C 1 by assumption, and so 
these are even smaller. Current observations constrain the curvature of the observable uni¬ 
verse to be \K/Hq\ < 0.01 [40]. Thus, the curvature terms in Eq. (5.8) are highly suppressed 
for any practical subhorizon calculation. 

We solve Eq. (5.8) as follows. We begin with linear solutions for the short scale variables 
6, V and 4>, i.e. neglecting Vp. We then correct short-wavelength scalar perturbations due to 
their coupling with long-wavelength modes. With this procedure, the solution, taking 5 for 
example, is 

5 = + (5,10) 

where is the linear solution, is the second-order correction due to mode coupling, 
and correspondingly for (j), u* and the velocity divergence. In this context it is simplest to 
work with the proper velocity divergence, 9 = Up^diV^ which takes into account the local 
scale factor. Further, we will switch from Tp to tp in the following. This is convenient, since 
eventually we would like to compare 6 with in the fiducial cosmology at fixed proper 
time tp. 

The growing mode for the potential is constant, i.e. (p = (piai- The linear (subhorizon) 
solution reads 

which obey equations for the fiducial cosmology with global expansion rate Ti and zero 
curvature Qm = 1- The continuity and Euler equations for the second order small-scale 
perturbations [i.e. those of the order given in Eq. (5.4)] then become 

j(2) + 0(2) = 

0(2) + 2i/0(2) + ^Lr25(2) = - 2afi\div{){djvi^^) 

-2 (Hp - H) 0(1) - 5 - H^) 5(1) , (5.12) 

where dots denote derivatives with respect to proper time tp. In the first line, we have used 
the fact that v\ is divergence-free due to its purely tidal nature. In Eq. (5.12), all coefficient 
functions should be evaluated at fixed proper time (which does not imply fixed conformal 
time as the scale factor is also modified). That is, Hp(tp) is the Hubble rate as function 
of time in the local modified cosmology, while H{tp) = 2/(3tp) is the same relation in the 
fiducial background. 

In order to solve for <5(2) and 0(2), we need initial conditions at early times Tp —?■ 0. In 
this paper, we will set the initial conditions to zero, (5(2)|ini = 0(2)|ini = 0, corresponding to 
no long-short mode coupling in the initial conditions. Ref. [10] showed that in single-field 
inflation in the attractor regime, there exists no local-type non-Gaussianity in the CFG frame, 
which implies that dj^) can only contain long-mode contributions of order In the scale- 

invariant case, they scale as {ks/'H)'^ [(^l/^s)^‘ h(ki)</)(k 5 )|ini]. These contributions, which 
are typically phrased in terms of “equilateral” and “orthogonal” bispectrum shapes, come 
from physical interactions during inflation, and are enhanced for example for an inflaton with 
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small sound speed. A calculation of the proper second order initial conditions including such 
terms is beyond the scope of this paper and will be considered in future work. Neglecting these 
terms is justified for “vanilla” single-field inflation with c* ~ 1, since these contributions are 
of order and thus suppressed by {l-L/ksY' (which eventually saturates 

at {'H/keq)'^) relative to the terms we are keeping. 

5.2 Initial conditions 

In order to solve for <5^^^ and 8^‘^\ we need initial conditions at early times Tp —>• 0. We are 
interested in deriving the late-time matter statistics for initial conditions provided by single¬ 
clock inflation. Thus, we need the small-scale density contrast in CFC in the presence of a 
long-wavelength mode, i.e. at second order, predicted by single-clock inflation. We obtain 
this by transforming the second-order result in the form of the bispectrum of the curvature 
perturbation TZ in comoving gauge as presented in [1] to CFC. This calculation has been 
presented in detail in [10], and we here summarize the results in order to make this section 
self-contained. 

In the notation of [1], but replacing Cthere with TZ in order to be consistent with our 
convention,' the long-wavelength metric perturbations in comoving gauge are given by 


500 — a^[—1 — 2A^i] 

9oi — ® Aj 

gij=a^[{l + 2TZ)6ij + jij], (5.13) 


where 'jij is tranverse-traceless and contains the tensor perturbations, which we will neglect 
in the following (see [21] for a CFC treatment of couplings between long-wavelength tensor 
modes and short-scale density modes [14, 45]). Note that we work to linear order here as 
we are only dealing with the long-wavelength mode, making use again of the advantages of 
the CFC formalism when dealing with a separation of scales. For the attractor solution of 
single-field inflation, the constraint equations in this gauge yield [1] 


_TZ _TZ' kl 

[aHY ■ 


(5.14) 


In this paper, we will ignore contributions from the initial conditions that scale as A:|, as 
they contain physical effects due to a nontrivial inflaton sound speed and thus depend on the 
particular inflationary model considered. We will discuss this momentarily. 

Neglecting A^i following the discussion above, the tetrad corresponding to the metric 
Eq. (5.13) is given by 


{eoY = a-\l,-N^) 

{eur = a-\l,[l-TZ]5^^). (5.15) 


Note that this implies that the proper time interval is dtp = a{T)dT. It is now clear that for a 
perturbation TZ constant in space and time (as both are suppressed by powers oikp), the only 
effect is to rescale the scale factor by a constant amount. The leading contributions in the 
small-fei, limit to the local Hubble rate are of order k'j^TZ. This suppression of both spatial and 
time derivatives makes physical sense: ii kp ^ Ti, then the local Hubble patch corresponds 

^Note that [10] followed [1] and used the symbol Z- 
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to a curved vacuum energy-dominated FLRW universe with curvature {kL/'H)‘^TZ. As the 
universe continues to inflate, with H ~ const, the curvature decays as and local Hubble 
patch rapidly asymptotes to flat de Sitter space (disregarding other perturbations that have 
been generated in the meantime). 

Given the Hp = H to the order in kp that we work in, and our freedom to rescale ap 
by a constant, we can simply choose ap = a and thus Tp = r. At lowest order in derivatives 
of the metric Eq. (5.13), the transformation to CFG at fixed t = Tp is given by = {ei)xp, 
that is 


X* = (1 - n)x^p . 


(5.16) 


Note that at this order in derivatives, the shift vector Ni drops out. Of course, we could have 
equally absorbed TZ into the definition of ap, yielding the same end result. 

The small-scale curvature perturbations r (to be distinguished from the long-wavelength 
perturbation TZ) then simply transform into GFG via 

rcFc('r,XF) = [l - TZxpdi] rgiobai(r, xj?). (5.17) 


This implies that the power spectrum of r transforms as 




1 + TZ 


din kgPx{ks) 
din ks 


[l + („^_l)7^]pglobal(^^), ^5 18 ) 


where — 1 is the slope of the dimensionless power spectrum k^P^, evaluated at kp, and 
Pr without superscript denotes the expectation value of the small-scale power spectrum 
(marginalized over long-wavelength modes). As shown in [1], the squeezed-limit bispectrum 
of curvature perturbations obeys the famous relation 

ks) = -{ns - l)PTi{kL)P.{ks) + o , (5.19) 

which is equivalent to a local modulation of the small-scale power spectrum by 


pglobal(^^) 


7 ^ 


[1 + (ns - 1)TZ] Pgi°^'^'(%). 


This modulation is exactly canceled by the transformation to GFG, i.e. 


pCFC 


(ks) 


7 ^ 


P,{ks) + 0{klTZ). 


(5.20) 


(5.21) 


That is, a local observer measuring the power spectrum of small-scale curvature perturbations 
during single-clock inflation cannot detect the presence of a long-wavelength mode up to order 
0{klTZ). 

Thus, for the calculation of small-scale perturbations in GFG frame, we have to set the 
initial conditions of the long-short mode coupling to zero, = 0. Gorrespondingly, = 

= 0. We emphasize again that for this we only needed to rely on a correct calculation 
of the leading long-short mode coupling (equivalently, the squeezed-limit bispectrum) in 
some global coordinates, along with the linear relation of the metric perturbation in these 
coordinates to the late-time density field. We do not need to rely on a second-order relation 
of P or r with a nonlinear definition of spatial curvature on some preferred spatial slice, as 
discussed in [41-43]. 
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Let us briefly discuss the corrections of order k\TZ oc that we neglect here. In the 
scale-invariant case, they scale as {ks[{kL/ks)‘^ These contributions, 

which are typically phrased in terms of “equilateral” and “orthogonal” bispectrum shapes, 
come from physical interactions during inflation, and are enhanced for example for an inflaton 
with small sound speed. A calculation of the proper second order initial conditions including 
such terms is beyond the scope of this paper and will be considered in future work. Neglecting 
these terms is justified for “vanilla” single-field inflation with ~ 1, since these contributions 
are of order <h(k/,)(/)(k 5 ) and thus suppressed by {'H/ks)'^ (which eventually saturates 

at {'H/keq)'^) relative to the terms we are keeping. 

5.3 Isotropic case: local expansion and cnrvatnre 

We first consider an isotropic long-wavelength perturbation, which implies = 0 and thus 
= 0. Restricting Eq. (5.12) to this case, taking a time derivative of the continuity 
equation and combining with the Euler equation, we derive one second order equation for 
5(2), 

5(^) + 2iL(ti7)5(2) - 5R'2(t^)5(2) = -2 [HpitF) - Hitp)] 5^^^ 

+ ^ [HUtF)n^{tF) - H^{tF)] (5.22) 

There are two approaches to relate the terms on the r.h.s. to the long-wavelength pertur¬ 
bations. One can either go back to the expressions in global coordinates [Eq. (3.15) and 
Eq. (4.14)]. This is detailed in App. B. Alternatively, one can use the fact that the scale fac¬ 
tor apitF) obeys the Eriedmann equation for a curved, matter-dominated universe, resulting 
in a Hubble rate given by 

Hliap) = H^fo + ^ko(^~f"] (5-23) 

where Hfo = HF{aF = 1) and = —Kf/Hpq. Note that is a constant, 

proportional to the comoving matter density. Evaluating this at early time, where Hf ^ 
H, > l,aF{tF) —> a{tF) [Eq. (3.18)], we obtain (evaluated at fixed If) 

H^F^ia% = H^a^ => HUtF)ni{tF) = H\tF)^^^ . (5.24) 

apitF) 

Now, is the physical matter density in the CEC region (up to a constant 3/(87rG)), 

while is proportional to the fiducial background matter density at the same proper time. 
Since the ratio between those matter densities is the density perturbation in synchronous 
gauge, we have 

= l + (5-25) 

apiTF) 

Note that this agrees with the r.h.s. of Eq. (3.19) after the difference between CEC time 
dtF = aF{TF)dTF and global time dt = a{T)dT is accounted for. We then obtain 

H‘j?{tF)^m{iF) - H‘^{tF) = H‘^{tF)^sc{tF) = ’ (5.26) 

where we have used Eq. (4.17) in the Einstein-de Sitter background (c)2<l> = S'H'^Asc/^). Note 
that this in fact holds to nonlinear order in Age [29]. 
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Now, going back to Eq. (5.24), we can evaluate this when Oi? = 1, which on the l.h.s. 
yields We can also evaluate it at a = 1, for which the r.h.s. yields Hq, which proves 

that = Hq. Thus, at linear order. 


Hpitp) - H{tp) = 


1 


2H{tp) 

= 2^iip) 




1 


T 2H{tp) 


[H^ap^ - Kpa-^ - H^a 


\ tp 


A 

^sc jj2 ® 

^0 J tp 


2 H{tp) ^2^ 
9 'H‘^{tp) 


With these results, Eq. (5.22) reduces to 

5(2) + _ 5 ^ 2 ^( 2 ) ^ J_g2^A ^(1) 

2 V 9 a2 J 


(5.27) 


(5.28) 


With the initial condition 5^‘^')\tp=o = 0 and = 0 imposed, we find the modified 

small-scale density field in the presence of the long-wavelength mode 

where all functions are evaluated at the same proper time. The second equality again uses 
Eq. (4.17). Eq. (5.29) agrees with solving spherical collapse in a homogeneous, but curved 
FLRW universe [35, 44], once the overdensity of that universe has been identified with the 
density perturbation in synchronous-comoving gauge. 


5.4 Anisotropic case: tidal field 

Following our discussion, we now set the isotropic contribution to zero, i.e. Hp ^ H and 
—?■ flm = 1 in Eq. (5.12), and keep only the terms involving Vp. We obtain the differential 
equations 

5(2) + 0(2) = _ a-i 

0(2) + 277 0(2) ^ 3^2^(2) ^ _ 2^-2 ^ . (5.30) 

An equivalent system of equations (within Newtonian gauge on small scales) was solved for 
EdS in [45]. Using their results [Eqs. (38), (39) and (45) there], we obtain 

= (5.31) 

where the anisotropic density tensor is defined as 

5.5 Velocity divergence 

We would also like to derive the velocity divergence 0 in CFG at fixed proper time tp. For 
this, notice that the continuity equation in Eq. (5.12), restricted to x* = 0 so that Vp = 0, 
reads 

e = -6. (5.33) 
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We can thus simply take the time derivative of Eq. (5.35) to obtain the second-order correction 
to 6 due to long-wavelength modes on the geodesic, 


0 ( 2 ) = H 


— A + 

21 ^ 



(5.34) 


The commonly used velocity divergence 9 defined with respect to comoving coordinates at 
this order is then simply 0(2) = a 0(2). 

The result for velocity divergence is useful because for pure potential flow (valid in 
the regime of perfect fluid) the velocity shear tensor can be derived from it. The velocity 
shear tensor is needed in order to derive the effects of redshift-space distortions on clustering 
statistics. 


5.6 Quasi-local Eulerian fields 

Combining Eq. (5.29) and Eq. (5.31), we write down the second-order small-scale density 
field quantifying the long-short mode coupling, strictly on the central geodesic, 

This describes the modified clustering induced by a long-wavelength gravitational potential 
that a local observer would measure along G relative to the average small-scale modes at fixed 
proper time. However, in order to measure clustering observables such as power spectra 
and bispectra, the observer also needs to survey how density varies across a finite region 
surrounding G. A computational trick is to choose the CEC central geodesic to pass through 
a given spacetime point under examination, and compute local physical quantities within 
that particular CEC. It is shown in Sec. 6 of [21] that measurements from neighboring CEC 
observers can be related by adding the local geodesic deviation. 

To be more specific, neighboring CEC observers are Lagrangian observers, since the 
coordinate system is tied to the fluid trajectory. Their measurements of the density field on 
each individual geodesic can be collected and mapped to a density field across a quasi-local 
Eulerian coordinate patch via 

Se = do — s* dido, (5.36) 


where the geodesic deviation vector s* is given by ® 



dr' 

apir') 



dT"aF{T") (didjhoo) , 


(5.37) 


where /i^g = —Kf^x'^x^, and the Lagrangian coordinate can be identified with x^ at the 
perturbative order we seek. For Einstein-de Sitter universe we find 

= ( 5 - 38 ) 

We will refer to the second term in Eq. (5.36) as the “displacement” contribution. The 
quasi-local Eulerian density field then reads 


dE = d^^^ + ^Asc <5(^) + . (5.39) 

®One does not have to distinguish between ap and a here because is first-order in the long mode. 
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In the same way, the velocity divergence (with respect to comoving coordinates) in quasi-local 
Eulerian frame, including the “displacement” term, reads 

Oe = - ‘^n Asc <5^^) - . (5.4o) 

Note that the explicit dependence on x* in the displacement terms will enter in the power 
spectrum as a derivative with respect to scale via the Fourier transform. Eqs. (5.39)-(5.40) 
are very similar to the expressions for the second-order density and velocity divergence in 
Eulerian Standard Perturbation Theory (SPT), which are restricted to subhorizon scales 
for all modes. Here we have generalized these results to a large-scale mode of arbitrary 
wavelength, albeit in the limit ki <C ks- The key differences to the SPT result are that first, 
the long-wavelength mode Age, is in synchronous-comoving gauge (since ks Ti, the 
gauge choice for the small-scale modes is not important). Second, the displacement term only 
enters as a relative displacement oc , which removes the gauge artifact in the unphysical 
absolute displacement present in the SPT expression. 

Eqs. (5.39)-(5.40) thus clearly show that intrinsic contributions to matter clustering due 
to long-short mode-coupling, which generate gravitational non-Gaussianity, are proportional 
to Age or Kfj, both of which scale as rather than the gravitational potential itself. 

This holds for long-wavelength modes on arbitrarily large scales. This has to be so, as the 
Equivalence Principle demands that constant or pure-gradient metric perturbations do not 
produce locally observable gravitational effects. Thus, any coupling to small-scale modes of 
the long-wavelength potential 4>(fci) itself has to be imprinted in the initial conditions. This 
does not happen in single-clock inflation. 

The expressions in Eqs. (5.39)-(5.40), representing locally measurable fields, do show 
some differences from the SPT results in terms of the numerical coefficients. This is due to 
two effects. First, within the local CFG patch, density fluctuations have been normalized 
with respect to the “local mean density” 

PFitp) = p{tF) [1 + Agc(tir)] , (5-41) 

which is modified by the long mode. Normalizing fluctuations with respect to the global p 
would result in an additional correction Agcd^^^ to the second-order density. Secondly, for 
both density and velocity gradient, the comoving coordinates (CFG) have been defined with 
respect to the local scale factor 

<iF{tF) = a,{tF) [1 — Asc(t_F)/3] . (5-4:2) 

Switching back to the global scale factor a amounts to providing an additional “displacement” 
term, which is equivalent to adding back the trace part to in the displacement s*. A 
detailed description of how to relate our results Eqs. (5.39)-(5.40) to SPT expressions can be 
found in App. D. Note that both of these additional corrections only involve Age oc 
and obey the Equivalence Principle. 


5.7 Squeezed matter bispectrum as measured by a local observer 

Following Eq. (5.39), the CFG observer measures a local power spectrum of short-scale density 
fluctuations (in the presence of long-wavelength perturbation). 


-P<5(k5)lA(ki) “ 




d\nPs{ks) 
din ks 
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P&{ks), 


(5.43) 
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where Ps{k) is the linear power spectrum of density fluctuations (in cN gauge, although this 
is not relevant at the scale ks)- Eq. (5.43) can then be used to calculate a “locally observ¬ 
able squeezed-limit matter bispectrum” as follows. Consider a class of observers distributed 
over large scales which measure the local mean density pp and small-scale power spectrum 
[l.h.s. of Eq. (5.43)] in their environment at fixed proper time. They then communicate this 
information to a distant observer in the future light cone of the entire class of observers. 
This distant observer in the future can then construct an estimate of the ensemble average of 
(A(ki)P5(k5)|^(kL))) which is equivalent to the squeezed matter bispectrum. The ensemble 
average of what he would measure is then 


(5(ki)5(k2)Asc(kL))' 




dlnPsjks) 
din ks 


P^{kL)Ps{ks)- 


(5.44) 


where k^ = (ki — k 2 ) /2 and psL = • k^ and we require |ki| ss |k 2 | ^ ks ^ kp. For an 

EdS universe, the linear matter power spectrum for the long mode in sc gauge is related to 
the primordial potential power spectrum via P^{kL) = {4:/97i^)kf^ P^{kL)- 

It is worth stressing that the scaling behaviour of the squeezed-limit bispectrum would 
be dramatically different for primordial non-Gaussianity of the local type [46] in the initial 
conditions. 


<hNG 


ini 




(5.45) 


where the non-Gaussian potential is a local function of a Gaussian potential 4>. In that 
case, the density field at second order is expressed through the first-order Gaussian fields via 


5 ( 2 ) 


ini 


2/^°£cI>5(i) 


(5.46) 


Instead of its second derivative, the gravitational potential is directly involved, and would 
lead to a contribution oc P^^^’^{kL) in Eq. (5.44), which leads to order unity departures 
from our result if /cl/H < (/nl)^^^- "^^e same effect gives rise to a scale-dependent bias 
for matter tracers on large scales [15]. From our analysis within the CFG formalism, we 
do not find an observable ® scale-dependent bias from general relativistic clustering, 

nor from a general relativistic specification of initial condition for single-clock inflation. Any 
locally observable effect of this type has to be imprinted in the initial conditions by a departure 
from single-clock inflation. This is different from the conclusions in Refs. [16, 17, 41, 43]. In 
combination with the results of [10] , this implies that in single-field inflation in the attractor 
regime, /^l— type contributions to the local matter bispectrum and halo bias vanish iden¬ 
tically; in particular, there is no f^f ~ (n* — 1) contribution to observables. Further, this 
constraint also applies to local-type non-Gaussianity of higher order, where higher powers of 
$ are added to Eq. (5.45), starting with 5(]^£<1>^. 

Of course, Eq. (5.44) does not correspond to how the matter bispectrum is measured 
in practice. In experiments, astrophysical tracers of the matter distribution are charted 
according to their observed redshift and their apparent position in the sky. Therefore, their 
number density, or the underlying matter density, appears distorted to a distant observer, 
since photons emitted from the source region are subject to change in the frequency and 
the direction of propagation as they travel from the source to the observer. Relevant effects 
include boundary effects such as redshift-space distortion and the Sachs-Wolfe effect, as well 
as line-of-sight integrated effects such as gravitational leasing, the Integrated Sachs-Wolfe 
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effect, and time delay. In a full general relativistic treatment, individual effects are gauge 
dependent, only the sum of them is gauge invariant. Detailed discussions of the projection 
effects are beyond the scope of this paper. We refer interested readers to a systematic study 
of the relativistic projection effects on the galaxy power spectrum [47] and on the squeezed 
galaxy bispectrum [10, 22, 48]. 

Projection effects can mimic the parameter in the apparent clustering. However, 
they arise purely as a result of photon propagation, and are physically unrelated to the nonlin¬ 
ear gravitational dynamics that occur locally (general relativistic or not). Projection effects 
in general are sensistive to the late-time expansion history as well as properties of source 
objects such as luminosity function and bias. Thus, by observing tracers with different bias 
and luminosity functions, it is in principle possible to measure the projection contributions 
separately and disentangle them from the local properties of the matter density field which 
we have derived here. 


6 Summary and conclusions 

In this work, we have proved that an inertial observer in a flat FLRW universe perturbed 
by a long-wavelength density perturbation (isotropic in the sense that didj^ = 6ijd‘^^/3) 
would locally experience a modihed FLRW expansion, with a renormalized scale factor ap, 
expansion rate Hp and an induced effective spatial curvature Kp. Those results have been 
rigorously derived by explicit construction of the CFC around observer’s geodesic. Since ap 
and Kp are uniquely determined as local observables (up to a multiplicative constant), there 
is no ambiguity in this construction. 

The necessary conditions for this “separate universe” picture are that the anisotropic 
stress vanishes and that sound horizons of all fluids are vanishing or negligibly small compared 
to the scale of the long-wavelength density fluctuation. We have presented expressions for 
Hp and Kp, at linear order in the long-wavelength perturbation amplitude and valid on 
arbitrarily large scales. The former includes a contribution from metric perturbations, in 
addition to peculiar velocity divergence. The latter is directly proportional to the Laplacian 
of the gauge-invariant curvature perturbation on comoving slices TZ. The CFC construction 
guarantees that these results are gauge-independent. Using this fact, the impact of long- 
wavelength density perturbations on small-scale dynamics can be numerically simulated by 
implementing the modified expansion rate ap and the modihed spatial curvature Kp [29]. 
This technique has a wide range of applications. For instance, the reponse of tracer abundance 
to a change in the “background” cosmology is related to tracer bias [35]. Another example 
involves power-spectrum response functions measured comparing separate universes with 
different FLRW parameters, which give the angle-averaged squeezed matter higher-point 
correlations [44]. 

Generalizing to anisotropic long-wavelength perturbations, the leading deviation ® from 
FLRW cosmology is a trace-free tidal contribution to the metric that is quadratic in the spa¬ 
tial CFC coordinates. The tidal tensor is exactly in the Newtonian form {didj — (5jj9^/3)<h 
without general relativistic corrections, if the scalar potential in conformal-Newtonian gauge 

is used. Since by construction CFC eliminates (essentially) all gauge freedom, the tidal 

^Corrections of higher order in the spatial departure from the observer’s trajectory would generate observ¬ 
able effects on small-scale dynamics that are suppressed by powers of the scale ratio These effects will 

have phenomenological relevance beyond the squeezed limit. 
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tensor is also independent of the gauge one chooses to parameterize long-wavelength pertur¬ 
bations. Implementation of the tidal tensor into separate universe simulations could be a 
potential techniqne to study tidal alignment of halo shape or galaxy shape. 

For an isotropic confignration of the long-wavelength pertnrbation, Birkhoff’s theorem 
guarantees that the “separate universe” can be generalized to all orders in the amplitude of 
the long mode along the observer’s geodesic (as proved in Sec. 2, for the tractable case of 
a compensated tophat). For the more general anisotropic case, the local expansion history 
apitp) is expected to violate the Friedmann equation at second order, as suggested by the 
presence of an extra term, quadratic in the anisotropic velocity shear, in the Raychaudhuri 
equation [49]. Despite the breakdown of the “separate universe” conjecture, our construc¬ 
tion of CFC will remain applicable. In that case, it would depict a locally uniformly and 
isotropically expanding spacetime, again up to tidal metric corrections. 

Another major result of this work is the squeezed-limit matter bispectrum measured by a 
local observer [Eq. (5.44)] for an EdS cosmology, which we derive using CEC without having 
to resort to the full second-order Einstein equations. While for simplicity the subhorizon limit 
is assumed for the short-scale modes, the full general relativistic dynamics of the long mode 
are acconnted for by the CEC formalism. Long-wavelength perturbations generate a squeezed 
bispectrum through their influence on the local gravitational clustering. This contribution 
takes the same form as in the familiar SPT result within Newtonian gravity, as long as the 
long-wavelength density contrast in sc gange is used. In other words, there are no explicit 
corrections dne to general relativity. This cnrious aspect at second order is analogous to the 
observation at linear order that the Newtonian Poisson equation can be extended to arbitrary 
scales if the sc-gauge density contrast and the cN-gauge potential are used [37] . In order to 
relate this bispectrum to observations from large-scale structure surveys, projection effects 
due to photon propagation through a perturbed background need to be included. While the 
ruler perturbations of [50] in principle provide all the necessary projection effects, we leave 
the detailed calculation of these to future work. 
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A Spatial gauge freedom 


In the presence of long-wavelength scalar perturbations, the space-space components of the 
CEC tidal metric (Eq. (3.23)) can be written as 


hfj = ^Kf {xiXj - - 6ij ^< 9 ^ 9 / 4 ' - ^ x^x^ 


(A.l) 


+ 


I {diidjdk + dijdidk + 2Sijdkdi - dkididj) 4 - - ^ {Sikdji Sudjk + dijdki) < 9 ^ 4 ' 
o y 


x^x^. 


In the first line, the first term is proportional to the local spatial curvature. The second 
term, proportional to 5ij, is equal to the same tidal potential in /iqo if 4^ = 4*. This is 
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the expected term satisfying the conformal Newtonian gauge condition hij = 6ij /iqo in the 
absence of anisotropic stress. Below we first show that the curvature term is equivalent to the 
stereographic form Eq. (3.27) for a curved FLRW universe. We then show that the second 
line can be eliminated by a re-parameterization of the spatial coordinates, and therefore has 
no physical consequences. 

We use the following spatial diffeomorphism: up to including C)[(x*)^], the most generic 
re-parameterization of spatial coordinates which leaves hoo and hot invariant reads [21] 

X* —> x'' + x^x'^xK (A.2) 

0 


The rank-4 tensor A^j^i is fully symmetric with respect to its last three indices. It is allowed 
to be a function of r, and the time-dependence only affects /iqo and hoi at C>[(x*)^]. The 
change to the line element is 

dx* = dx* -|- ^A^jki x^x’' dx\ (A.3) 

which changes hij by 


^ij ^ ^ij T 2 T ^ ^ • 


(A.4) 


A.l The effective curvature 


We have demonstrated in Sec. 4 that long-wavelength scalar perturbations must induce an 
effective curvature Kp in the local-coordinate description. Here we see that the effective 
curvature is residing in the first term of Eq. (A.l). 

Now in Eq. (A.2) we choose 

A^jkl = {^^j^kl + + ^^i^jk) , (A.5) 

which generates radial reparameterization Eq. (3.26). This changes that metric term by 

^Kp [xiXj - 6ijr'^) —^ -^Kpdijr'^, (A.6) 

and exactly accounts for the lowest order curvature contribution of the stereographic param¬ 
eterization Eq. (3.27). 

It can be concluded that as long as Kp is conserved, the only local isotropic effect of 
long-wavelength scalar perturbations, at ©[(x)^)^], is an effective FLRW curvature, with no 
further general relativistic correction. 


A.2 Conformal Newtonian coordinates for CFC 

What about the anisotropic tidal terms? Since we only have rank-2 tensors didj'^, and 
5ij at our disposal, the most generic Aijki that has the right symmetry can be written as 

Aijkl — ^ i^dijSkl “I” hikdjl -\- hnSjk] d ^ b (^Sijdk^l T ^nd^dj “t" hikdldj') T 

+c (hkididj + djididk -h Sjkdidi) T, (A.7) 

where a, b and c are arbitrary numerical coefficients. By matching the coefficients of various 
terms, we find a unique solution in order for the change in Eq. (A.4) to exactly compensate 
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the second line of Eq. (A.l), a = —1/9, b = 2/3, c = —1/3. With a successive change of 
coordinates in Eq. (3.26), we arrive at the following simple form for the CFC metric 


ds'^ 


apirp) 


- 1 + dkdi<^ 




dTp 








Sijdx^dx^ 
{1+Kpr‘^/Af _ ' 


(A.8) 


We therefore have put the CEC metric into the conformal Newtonian form, which reduces 
exactly to the Newtonian coordinates on sub-horizon scales. It is manifest in this conformal 
Newtonian frame that all the anisotropic effects of the long-wavelength perturbation are 
parameterized by a simple tidal potential term. 


B CFC expansion rate and matter density 

Applying Eq. (3.15) and Eq. (4.14) to EdS universe with growing perturbations, we have 
Hp{tp) - H{t) = -h(^^ + 

HUtp) - H\t) = ^ + HUtp) - H\t) = -2HH + (B.2) 

The global time t is understood to refer to the same spaeetime point which has proper time 
tp in CFC. Therefore, t and tp are numerically different, which can only be neglected on the 
right hand sides of the equations which explicitly contain first-order perturbation variables. 
We furthermore need to account for modification of the proper time due to long-wavelength 
perturbations, and derive H[t) — H[tp) = —H At, with At = tp — t = 2^/3H, and similarly 
for H‘^{t) — H‘^{tp). We eventually find 

2JJ 

Hp{tp) - H{tp) = (B.3) 

Hp{tp) ^^{tp) - H‘^{tp) = (B-4) 

Compared with Eqs. (B.1)-(B.2), we see that all terms proportional to the Newtonian poten¬ 
tial itself cancel out; locally observable effects must depend on second spatial derivatives 


C The 2nd Friedmann equation 


For completeness, we should also verify that the second Friedmann equation holds, which 
can be done in an analogous way. It follows from Eq. (3.20) that the acceleration is given by. 


1 d'Hp 1 uap 


ai drp 


ap dt'j^ 


1 

a dP 


I 


+ 2'H'4> + -HT' + + nd-V) 

o 


(C.l) 


where dtp = ap{Tp)dTp is the CFC proper time. We now use the trace (isotropic) part of 
the space-space Einstein equation (note that there is in general a total pressure perturbation 
on the right-hand side), 


2T" + n (4T' + 2$') + 2 {2n' + n^)^ = 8pGa^ {Vp - V) , (C.2) 
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We then have 


n^' + 2%'^ + n^' + ^” -^{d-V + nd- v) 

3 

= \ [02^ - 3n (^' + n<^>)] + 47rGa2 {Vf - P) + \d‘^ (4> - ^') 

o o 

= \{PF — p) + 3 {Vf — Py\ + (‘h — 4'), (C.3) 

via the time-time Einstein equation. Again, the second term vanishes if there is no anisotropic 
stress. Inserting this into Eq. (C.l), we immediately obtain the form of the second Friedmann 
equation, 


1 d?aF 
ttF dPp 


AfG 


{pF + 3Vf) , 


(C.4) 


where again pF and Vf are respectively the total energy density and pressure at the observer’s 
location. Contrary to the discussion of the first Friedmann equation, no assumption on the 
cosmic fluids is needed for the second Friedmann equation to hold. 


D Relating CFC result to SPT result 


How does Eq. (5.39) compare to the usual result of Standard Perturbation Theory (SPT) at 
second order? First of all, the density contrast measured in CFC is computed with reference 
to the “local homogeneous density”, i.e. pF = /o(l + Age), while in SPT one makes reference 
to the global homogeneous density p. This amounts to adding back a contribution Age <5^^^. 
Secondly, the “displacement” term in Eq. (5.39) is proportional to the trace-free Kf^ This is 
because the trace part, corresponding to an isotropic convergent or divergent flow, has been 
absorbed into the local scale factor qf- To undo this renormalization, one adds back a term 
of this isotropic flow (to find aF{tF)/a{tF) — 1 one uses Eq. (3.19) and then accounts for 

At = If — t), 


/ ajtF) _ A 

\aF{tF) J 




1 

3 




Implementing these two corrections, we obtain the second-order density constrast 


,5(2) = 


didj 
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^ Age c5(^) + ^ 
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d^d^ 

d2 


sW 


(D.l) 


(D.2) 


Finally, the spatial origin is still chosen at the center of the local patch, and is therefore 
moving with it. One can undo the bulk motion of the local patch by including a homogeneous 
“displacement” (which leaves correlation functions invariant), and then adjust to whatever 
new spatial origin that is desired. 


(^Age^ x^di6^^'> ^ [1 + x^dj] (^Age^ di5^^^ ^ 

We obtain 


= (A A. 
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d2 


d^d^ 
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(D.4) 
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These tensorial structures and numerical coefficients agree completely with the familiar F 2 
kernel in second-order SPT computed for EdS universe [51], in the limit that one wave 
number is much larger than the other, except for an overall factor of two, which accounts 
for swapping the long-wavelength and short-wavelength modes. We stress that the F 2 kernel 
is derived from Newtonian nonlinear clustering on subhorizon scales, while our result is still 
valid for (super-)horizon scale ki. The Newtonian form is, however, preserved when the 
long-wavelength density contrast in sc gauge is used. 

Performing a similar check on the velocity divergence Eq. (5.40), we find agreement 
with the G 2 kernel of SPT. 
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